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Abstract 

In this paper the entanglement and quantum phase transition of the anisotropic s = 1/2 XY 
model are studied by using the quantum renormalization group method. By solving the renormal- 
ization equations, we get the trivial fixed point and the untrivial fixed point which correspond to 
the phase of the system and the critical point, respectively. Then the concurrence between two 
blocks are calculated and it is found that when the number of the iterations of the renormalzia- 
tion trends infinity, the concurrence develops two staturated values which are associated with two 
different phases, i.e., Ising-like and spin-fluid phases. We also investigate the first derivative of the 
concurrence, and find that there exists non-analytic behaviors at the quantum critical point, which 
directly associate with the divergence of the correlation length. Further insight, the scaling behav- 
iors of the system are analyzed, it is shown that how the maximum value of the first derivative of 
the concurrence reaches the infinity and how the critical point is touched as the size of the system 
becomes large. 
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I. INTRODUCTION 



The quantum entanglement, as one of the most intriguing features of quantum theory, has 



attracted much attention because its non-classical correlation can 
resource in quantum communication and information processing 
tion between the quantum entanglement and quantum correlation 



De regarded as an essential 
ll. In view of the connec- 
the entanglement in 



some many-body systems has been widely investigated . The entanglement, which is 
the main difference between quantum and classical systems, also plays an important role 
in the quantum phase transition (QPT). The QPT is induced by the change of an external 
parameter or coupling constant 6j. This change occurs at absolute zero temperature where 
the quantum fluctuations paly the dominant role and all the thermal fluctuations get frozen. 

The entanglement between two nearest-neighbor spins in one-dimensional XY system 
with the transverse magnetic field was studied in Ref. The relation between the QPT 
and entanglement are studied in this paper and it is found that the system exists QPT and 
the scaling behavior in the vicinity of the critical point. After that, a great deal of efforts 
have been devoted to investigating the entanglement and QPT in some spin systems js-ll]. 
The XY model is exactly solved by the Jordan- Wigner transform in Ref. js], in which the 
Ising model, as a special case of the XY model, exhibits QPT and a maximum value for 
the next-nearest-neighbor entanglement at the critical point. The entanglement properties 
and the spin squeezing of the ground state of mutually interacting spins 1/2 in a transverse 
magnetic field are analyzed and the system shows a cusplike singularity at the critical point 
in the thermodynamical limit \lX]- To further discuss the QPT and entanglement, the 
renormalization group (RG) method is introduced, such as real space RG, Monte Carlo RG, 



and density-matrix RG 
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15| . The pairwise entanglement of the system is also discussed 



16 



17|. Very recently, the 



by means of quantum renormalization group (QRG) method 
spin— 1/2 Ising and Heisenberg models are studied by using the same method by a group of 
Iran and found that the systems exist QPT It is also shown that the nonanalytic 

behavior of the entanglement and the scaling behaviors closing to the quantum critical point 
are obtained. 

The XY model is firstly discussed in Ref. |22| , and then has been widely investigated re- 



cently 



23N25|. In this paper, we study the quantum entanglement and QPT in the spin— 1/2 



XY model by the method of QRG. It is found that the system exists QPT between the spin- 
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fluid and Ising-like phases. Further insight, the nonanalytic behavior of the entanglement 
and the the scaling behavior of the system are gotten. The organization of this paper is as 
follows. In Sec. [Tllwe briefly introduce the RG method and obtain the fixed points of the 
XY model. The block-block entanglement is obtained in Sec. Illll and in Sec. IIVI we discuss 
the non-analytic and the scaling behaviors of the entanglement. Sec. |V]is devoted to the 
conclusions. 



II. RENORMALIZATION OF THE XY MODEL 

Firstly, the QRG is introduced. Eliminating the degrees of freedom of the system followed 
by an iteration is the main idea of QRG. The aim of the iteration is that reduces the number 
of variables step by step until a more manageable situation is reached. For this idea, the 
kadanoff's block approach is implemented in this paper, because it is not only well suited 
;o perform analytical calculations but also easy to be extended to the higher dimensions 
We consider three sites as a block, the kadanoff's block approach is given in Fig. 1. 
In this way, the effective Hamiltonian (if^-^-^) can be gotten which has structural similarities 
with the original Hamiltonian {H). The original Hilbert space is also replaced by a reduced 



Hilbert space which acts on the renormalized subspace 271. |28[|. 



The Hamiltonian of the XY model on a periodic chain with N sites can be written as 

T ^ 

H(j^^) = jYl [(1 + 7) ^r^r+i + (1 - 7) ^r^i] , (i) 



where J is the exchange coupling constant, 7 is the anisotropy parameter, and a" (a = x, y) 
are Pauli matrices. The XY model can be encompassed another two well-know spin models, 
i.e., the Ising model for 7 = 1 and the XX model for 7 = 0. For < 7 < 1, it belongs to 
the Ising universality class [22]. 

By using the Kadanoff's block approach, Eq. ^ can be written as 

H = H^ + H^^, (2) 

where is the block Hamiltonian, H^^ is the interblock Hamiltonian. The specific forms 
of and if are 

Ar/3 

H^ = Y.ht (3) 

L 
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N/3 

= I E + ^) ^1,3^2+1,1 + (1 - 7) 4,3^+1,1] , (4) 

L 

where 

= ^ [(1 + 7) Ki^I,2 + ^2,2^2,3) + (1 - 7) + 4,ai,3)] , (5) 

which is the Lth block Hamiltonian. 



In the terms of matrix product states [29|, the Lth block Hamiltonian can be exactly diag- 
onalized and solved. We can obtain four distinct eigenvalues which are doubly-degeneracy. 
Defining |'|') and \D as the eigenstates of a^, both degenerate ground states are given as 
follows: 

i$o) = ^ i — - (-VT+7\m) + V2 int) - i;n) + v^7 m)) , (6) 
$;) = ^ i — - (- v^7 ittt) + im) - i;n) + yr+y lut)) • (7) 

The energy corresponding the ground states is 

J 



^o = -^Vl + y. (8) 

For eliminating the higher energy of the system and retaining the lower, the projection 
operator Tq is composed by the lowest energy eigenstates of the system. The relation be- 
tween the original Hamiltonian and the effective Hamiltonian can be given by the projection 



operator 30|, i.e., if*^-^-^ = TqHTq, where Tq is the hermitian operator of Tq. In the ef- 
fective Hamiltonian, we only consider the first order correction in the perturbation theory. 
The effective Hamiltonian is 

if-// = H^ff + Hl^^ = T+H^To + T^H^^Tq. (9) 

The projection operator Tq can be searched in a factorized form 

To = ^fiTo^ (10) 

where Tq is the Lth block, which is defined as 

T^ = \^)l{%\ + W)l{%\- (11) 

The Ifl^)^ and in the Eq. f ill I) are the renamed states of Lth block which can be seen 
as a new spin— 1/2. The renormalization of Pauli matrices are given by 

T^oXl^o^ = ^r^L (^ = 1,2,3; a = x,2/), (12) 



where 



1 + 7 (1 + 7)' 



VW+l^' 2(1 + 7 



= /2(l77-)' '^' = -WT7)- ^''^ 



Then, the effective Hamiltonian of the renormahzed chain can be gotten, 

, 7V/3 

= T ^ [(i + V) ^^^^+1 + (i - V) ' (14) 

i 

where 

In the above equations, not only do we get the nontrivial fixed point 7 = 0, but also obtain 
the trivial fixed point 7 = 1 which specifies the critical point of the system by solving 7 = 7'. 
When 7 — )■ 0, the model falls into the universality class of XX model corresponding to a 
spin-fluid phase; while for 7 — )■ 1, the system is in Ising-like phase. From above analyzing, 
it is found that there is a phase boundary that separates the spin-fluid phase 7 = 0, from 
the Ising-like phase < 7 < 1. 

III. ENTANGLEMENT ANALYSIS 

We investigate the ground-state entanglement between two blocks of the XY chain by 



the concept of concurrence 



31 



32| and demonstrate how the concurrence varies as the size 
of the block becomes large. We consider one of the degenerate ground states to define the 
pure state density matrix. Thus, the pure state density matrix is defined by 

P=|$o)($o|, (16) 

where |$o) has been introduced in Eq. ([6]). The results is same if we consider the Eq. 
(I7j) to construct the density matrix. Because the concurrence is one measure of pairwise 
entanglement, we must trace over the degrees of freedom of one site in the block. Without 
loss of generality, we trace over the site 2. The reduced density matrix for the sites 1 and 3 
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can be obtained as, 



Pl3 



4(72 + 1) 



1 



/ 2 27 \ 

72 + 1 72 + 1 

72 + 1 72 + 1 

\27 272/ 



(17) 



Ci3 denotes the concurrence of the sites 1 and 3 which is defined as 



C 



13 — 



Max<^ A 



1/2 .1/2 .1/2 ,1/2 ^ 



} 



(18) 



where the (/c = 1, 2, 3, 4) are the eigenvalues oi R = pisfhs [where pis = 
{af (T3) {af erf)] in ascending order. The eigenvalues of R can be exactly solved 
which are 



We substitute the above equations into Eq. flTK]) . then C13 can be gotten which is the 
function of 7. The renormalization of 7 defines the evolution of concurrence on the increasing 
of the size of the system. For better to discuss the concurrence, we set C13 as a function of 
g, where 



Ci3 versus different QRG iterations is plotted in Fig. 2. 

For different QRG steps, the plots of C13 versus different QRG steps cross each other 
at the critical point. In the thermodynamic limit, the concurrence develops two saturated 
values which is nonzero for Qc = ^ and zero for < g < 1 and g > 1. At the point gc = 1, the 
system exists quantum correlation because C13 is a nonzero constant. The infinite chain can 
be effectively described by a three-site model with the renormalized coupling constants. In 
this case, the quantum fluctuations play an important role and destroy any long-range order 
of the system. Nonzero C13 verifles that the system is entangled where the ground state 
is characterized by a gapless excitation and algebraic decay of spin correlations. While for 
< g < 1 and g > 1, the model is gap with magnetic long-range order, the nontrivial points 
correspond to two Ising phases ordered in y direction {g — )• 0) and x direction {g 00), 
respectively. It can be seen that the Neel phase is the dominant phase of the system. 



Ai — A2 — 0, A3 




(19) 



(? = (l+7)/(l-7). 



(20) 
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IV. NON-ANALYTIC AND SCALING BEHAVIOR 



The system can occur QPT because of the non-analyticity property of the block-block 
entanglement. The non-analyticity behavior can be also accompanied by a scaling behavior 
because of the diverging of the correlation length. In this section, we show the behaviors of 
the concurrence and the QPT in the XY model. By using of QRG method, a large system 
(A^ = 3"^^) can be effectively described by three-site block with the renormalized coupling 
constants after the nth iteration of RG. Thus, the entanglement between two renormalized 
sites represents the entanglement between two blocks of the system each containing 
sites. 

The first derivative of concurrence is analyzed and shows a singular behavior when the 
size of the system is infinity. All these dates have been shown in Fig. 3. It is easy to see that 
the first derivative of the concurrence is discontinuous at Qc = 1, while C13 is continuous. 



This indicates that the QPT of the system is the second-order QPT [33|. The scaling 
behavior of the maximum oi y = dCi^/dg versus is plotted in Fig. 4 which shows a linear 
behavior of In(ymax) versus In(A^). In order to show the intuitionistic scaling behavior, the 
numerical results are given. The position of maximum {gm) of y touches the critical point 
as the size of the system increases. This is plotted in the Fig. 5. It is found that there are 
the relation of Qm = ~ with 6 = 0.98. The exponent 6 which is called entanglement 
exponent is directly related to the correlation length exponent closing to the critical point. 
The singular behavior of the concurrence and the scaling behavior of the system depend on 
the entanglement exponent, which is the reciprocal of correlation length exponent, i.e., 6 
= The emerging singularity connects to the universality class of the model. As the 

critical point is approached in the limit of large size, the correlation length covers the whole 
system. That is to say, the RG implementation of entanglement truly captures the critical 
behavior of the XY model in the vicinity of the critical point. 



V. CONCLUSIONS 



In this paper, we discuss the entanglement and QPT in the anisotropic s = 1/2 XY 
model by the quantum renormalization group method. The concurrence as one measure of 
the quantum correlation is investigated. The critical behavior of the system is obtained by 
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the renormalization of the lattice. As the number of RG iterations reaches the infinity, the 
system occurs QPT between the spin-fluid and the Ising-like phases which correspond to two 
different fixed values of the concurrence at the critical point and both sides of it respectively. 
The diverging behavior of the first derivative of the concurrence is accompanied by the 
scaling behavior in the vicinity of the critical point. Further insight, the scaling behavior is 
investigated and characterizes how the critical point of the model is touched as the size of 
system increases. In the thermodynamic limit, the non-analytic behavior of entanglement 
is correlated with the diverging of the correlation length at the critical point. 
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Figure captions: 

Fig. 1. The kadanoff 's block renormalization group approach of a chain where we consider 
three sites as a block. 

Fig. 2. Representation of the evolution of the concurrence versus g in terms of QRG 
iterations. 

Fig. 3. First derivative of concurrence and its manifestation toward diverging as the 
number of QRG iterations increases (Fig. 2). 

Fig. 4. The logarithm of the absolute value of maximum, InddCis/ dg\giaax\) , versus the 
logarithm of chain size, ln(A'"), which is linear and shows a scaling behavior. Each point 
corresponds to the maximum value of a single plot of Fig. 3. 

Fig. 5. The scaling behavior of gimax in terms of system size (N) where gimax is the position 
of maximum in Fig. 3. 
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